Relativistic Quark Model Calculation of the l\, l 2 Coefficients of the Chiral Lagrangian 
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We briefly report on a relativistic quark model scheme to calculate the 0(P 4 ) pion-pion vertex 
in the planar approximation and in the chiral limit. The calculation is reduced to the solution of 
simple integral equations (Bethe-Salpeter like) by an effective use of chiral Ward Identities. Specific 
model computations are provided. 
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The interest of meson scattering in the quark model has recently revived due to the exotic 1 
BNL at the energies of 1.4 and 1.6 GeV . The hybrid @ or meson molecule nature of these states is being debated. 
The standard tools for meson scattering are the chiral lagrangian |^) and its unitarized extensions ||, which cannot 
access that intermediate energy region, let alone interpret the Fock space nature of the resonances. Still they provide 
the benchmark for model approaches to the underlying dynamics which encompass chiral symmetry. In recent works 
||, ||] the use of a novel chiral Ward Identity (|3| below) allowed the reduction of the planar quark model approximation 
to pion scattering, 
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to reproduce the 0(P 2 ), 0(m%) chiral amplitude. Here we continue this effort by establishing contact with the 
0(p 4 ), 0(m°) chiral lagrangian. This work is related to microscopic studies of it — it scattering with the Nambu and 

"]L with the equal time quark model [Q, ^] , and with a direct 
Cfl . While the 0(P 2 ) amplitude is model independent, the 
0(p A ) provides a model test. Here we exactly and exp icitly determine the Feynman diagrams that govern the 0(p 4 ) 
amplitude, then employ two simple models for a test evaluation. 

In the standard notation, the pion scattering amplitude arising from its 0{P A ) contact terms alone reads 



Jona-Lasinio model |?j , with the bosonization technique 
numerical evaluation with a Schwinger-Dyson model [ I| 
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The Bethe-Salpeter pion vertex \-n in (0) above, function of two momenta (incoming meson P and quark k + P/2) can 
be conveniently expanded in an external momentum power series. But first it is useful to extract from it a dressed 
vertex Ta satisfying interesting Ward Identities connecting it to the fermion dressed propagator and the planar ladder, 



x{P>k) = 



-iT A {P,k) + A(P,k) 
V2U 



(3) 



T A (P,k + P/2) = S- 1 (k-P/2) l5 + l5 S-\k + P/2) 
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The use of these identities, as well as the fact that no zero-momentum pole appears in any vector or scalar ladder 
given the symmetry breaking pattern of the theory, allows us to reduce the momentum four term of eq. (Q) to three 
different classes of diagrams, parameterized as, 
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All external momentum and color dependence is explicit. To calculate the ladders which remain in eq. (Q), we 
write integral equations for them, analogous to the Bethe-Salpeter equation for a meson and with up to ten indepen- 
dent functions of q 2 . Matching this Feynman amplitude with the physical external momenta and multiplying each 
permutation by the adequate isospin factors, we obtain for the Vs defined in eq. Q), 



h = — (-2d 1 +d 2 -6d 3 + 2d 4 + 5d 5 + d 6 ) 
3 

h = - — (d 2 -2d 3 + 2d 4 + d 5 -d 6 ). 
Id 
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The necessary Dirac gamma traces are computed with FORM || . We perform model evaluations in Euclidean space 
with an interaction rung given by V = —g 2 "f fl iK(k — q)^ which is similar to the One Gluon Exchange interaction 
in the Feynman Gauge. We have used two different infrared strong kernels, both UV finite to avoid renormalization 
complications. 

The first is a simple Gaussian kernel (whose Euclidean angular integrals are simply Bessel functions) K (q) = 
exp(—q 2 /A 2 ) with parameters g = 3, A = 1.15 GeV (model 1) which leads to a constituent quark mass of 116 MeV, a 
condensate —(170 MeV) 3 , a pion decay constant f n = 87 MeV and gives l\ = —0.11, l 2 — 0.14; while the parameters 
(model 2) g = 1.7, A = 2 lead to a mass of 93 MeV, condensate -(229MeV) 3 , and / w = 86 MeV, h = -0.06, 
l 2 = 0.12. The second is a rational kernel (whose Euclidean angular integrals contain only polynomials and square 
roots) K(q) = [l/(q 2 -\ 2 )-l/(q 2 -A 2 )} with g = 10.3, A = 0.8 GeV, A = 0.85 GeV (model 3), we obtain a constituent 
quark mass of 216 MeV, condensate —(156 MeV) 3 , pion decay constant 92 MeV , and Z 4 = —0.065, l 2 — —0.12. 

The computed l\ and l 2 have the correct signs and ratio when compared with the phenomenological analysis, 
however they are large. Investigations are underway to test a wider class of quark-antiquark interactions. 
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FIG. 1: Mass gap function solution of the Schwinger-Dyson equation for model 1 (in GeV). 
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FIG. 2: Solution of the homogeneous Bethe-Salpeter equation. From top to bottom, F0[o], E2[q], G0[q] (corresponding to the 
coefficients of an external momentum expansion of the llewelyn- Smith wavefunctions). 




FIG. 3: Solution of the homogeneous equation for two A's with ladder, needed to calculate dl,d2,d3,d4. (Coefficients of external 
momentum expansion). 
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FIG. 4: Solution of the homogeneous equation for P with ladder, needed to calculate d5,d6. (Coefficients of external momentum 
expansion). 



